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ABSTRACT 

We discuss possible distortions of the ionization history of the Universe in the 
model with small scale baryonic clouds. The corresponding scales of the clouds are 
much smaller than the typical galactic mass scales. These clouds are considered in a 
framework of the cosmological model with the isocurvature and adiabatic perturba¬ 
tions. In this model the baryonic clouds do not influence on the cosmic microwave 
background anisotropy formation directly as an additional sources of perturbations, 
but due to change of the kinetics of the hydrogen recombination . We also study the 
corresponding distortions of the anisotropy and polarization power spectra in connec¬ 
tion with the launched MAP and future PLANCK missions. 
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1 INTRODUCTION. 

One of the most important problems of the modern cosmol¬ 
ogy is the determination of the density and spatial distribu¬ 
tion of the baryonic fraction of the matter. 

There are several sources of information about Ulbh 2 = 
pb/Per parameter, where pt> and p cr are the present values 
of the baryonic and critical densities and h is the Hubble 
constant normalized to 100 kms _1 Mpc _1 . Firstly, the bary¬ 
onic fraction of the matter manifests itself in the well known 
mass -luminosity relation for galaxies and cluster of galax¬ 
ies which leads to the following value of the Ulbh 2 param¬ 
eter: Clbh 2 ~ 0.028lo oog (see for the review by Freedman 
et al. 2001). Another one comes from the confrontation of 
the Standard Big Bang Nucleosynthesis (SBBN) theory and 
observational data (see for the review by Fukugita, Hogan 
& Peebles 1998). The corresponding value of the baryonic 
density from this method is Vtyh 2 = 0.019 ± 0.001. An addi¬ 
tional empirical relation between baryonic and dark matter 
fractions Fb,m = — 0.1 at h = 0.65 comes from X- 

ray data on clusters of galaxies (Carlberg et al. 1996; Ettori 
& Fabian 1999). For the most popular ACDM cosmologi¬ 
cal model with ~ 0.3 and ~ 0.7, the corresponding 
value of the Ulbh 2 parameter is ~ 0.02 in agreement with the 
SBBN predictions. 

An independent important information about the bary¬ 
onic fraction of the matter in the Universe comes from 
the recent CMB experiments such as BOOMERANG (de 


Bernardis et al. 2000) and MAXIMA-1 (Hanany et al. 2000). 
Fitting the CMB anisotropy power spectrum to the above 
mentioned observational data (Tegmark & Zaldarriaga 2000; 
White et al. 2000 and Lesgourgues & Peloso 2000) indicates 
that a baryon fraction parameter should be significantly 
larger than the SBBN expected value, namely, Qbh 2 ~ 0.03. 
However, recently Bond & Critteden (2001) show that new 
BOOMERANG, MAXIMA-1 and DASI data do not contra¬ 
dict to Qbh 2 = 0.022 ± 0.004. 

It is worth noting that the above mentioned methods 
of the baryonic fraction density estimation from the CMB 
and SBBN predictions are based on the simple idea that 
the distribution of matter (including dark matter particles 
and baryons) is practically homogeneous for all scales, ex¬ 
cept some fluctuations leading to the galaxy and large-scale 
structure formation. Typically, they are assumed to be adi¬ 
abatic one. One can ask, how sensitive are the CMB data 
themselves to the presence of the small- scale baryonic (non¬ 
linear) clouds before cosmological recombination and how 
can they transform the standard schemes of the cosmolog¬ 
ical parameter extraction from the CMB data? Definitely, 
this possibility is related to the isocurvature perturbations 
of the composite fluid which contains baryons, CDM parti¬ 
cles, photons and neutrinos at very high redshift z 2> 10 3 . 
There are a lot of modes of perturbations in the compos¬ 
ite fluid, which is discussed by Riazuelo & Langlos (2000), 
Bartolo, Matarrese & Rioto (2001), Polarski & Starobinsky 
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(1994), Abramo & Finelli (2001), Bucher, Moodley & Turok 
(2000) and others. The general idea about classification of 
modes of perturbations is based on a very simple definition 
of the isocurvature modes. They do not perturb the gravi¬ 
tational potential. This means that the fluctuations of the 
total matter density ptot are zero (see Burns 2001), 

N „„ 

C/Tl 

(5 ptot — ^ ^ pi$i + 4p-,.(l T Ris'y') j-, — 0, (1) 

i =0 

where pi denotes the density of each massive species in¬ 
cluding baryons and different kinds of the CDM particles, 
Si = Spi/pi is the density contrast for each massive com¬ 
ponent, R V1 is the density ratio between neutrinos p v and 
black body radiation p 7 , and ST/T is the CMB temperature 
perturbations. 

We would like to point out that in the definition of 
the isocurvature perturbations in Eq. (|l|) one can find some 
peculiar mode (or modes) which compensates the baryonic 
perturbations potential, i.e., it corresponds to the condition 
phSb = —p x 8x for some x component of the CDM particles 
mixture. We call below this mode as a compensate isocur¬ 
vature mode (CIM) for the ^-component of the dark matter 
particles. If several components of the dark matter parti¬ 
cles take part in the CIM formation, we will continue to call 
them as x- component.^] 

In principle, for the CIM perturbations it is possible to 
assume that the amplitudes 8 X and 5 b are less than unity or 
|(5a,| ~ 1 and |<56| ~ — {px/ Pb)\8 x \ 1- One of the most inter¬ 
esting cases corresponds to the model with 8 X — —1, which 
means that some patches of the cosmological matter do not 
contain the CDM ^-particles, but at the same patches there 
are non-linear clouds of the baryonic matter which compen¬ 
sate the perturbations of the gravitational potential. Below 
we will assume that a typical mass scale of the CIM pertur¬ 
bations is smaller than the typical galactic mass scale. This 
means that CIM perturbations do not influence on the CMB 
anisotropy formation as the additional sources of perturba¬ 
tions, but they can transform the kinetics of the hydrogen 
recombination. This leads to the transformation of the cor¬ 
responding Ci power spectrum of the CMB for the adiabatic 
fluctuations at the scales above a few Mpc. 

It is necessary to note that the idea about non- 
homogeneous distribution of the baryonic matter at small 
scales is not new. The importance of the entropic pertur¬ 
bations in the history of the cosmological expansion was ad 
hoc demonstrated by Doroshkevich, Zel’dovich and Novikov 
(1967) and Peebles (1967,1994) and recently was generalized 
taking into account multi-species structure of the cosmolog¬ 
ical plasma by Gnedin & Ostriker (1992), Hogan (1993) & 
Loeb (1993), Peebles & Juszkiewich (1998). The possible 
inhomogeneities of the baryon fraction distribution in the 
epoch of the nucleosynthesis (Inhomogeneous Big Bang Nu¬ 
cleosynthesis -IBBS) was widely discussed in the literature 
(see for the review by Jedamzik & Rehm 2001) in connection 
with quark-hadron phase transition. But the typical scales 

* As one can see this mode corresponds to 5T/T = 0. This means 
that the CIM are equivalent to an isotemperature perturbations. 
Note, that exactly the same mode was described by Abramo & 
Finelli (2001), but for compensation between quintessence scalar 
field perturbations and some kind of the CDM-particles. 


of such kind of peculiarities are extremely small compared 
to the typical mass scale M ~ 10 5 — 10 6 Mq for the isocur¬ 
vature perturbations. Other events or processes have been 
suggested as possible sources of the isocurvature perturba¬ 
tions partly connected with the baryon re-distribution in the 
space. For example, cosmic strings and corresponding cur¬ 
rents and magnetic fields could generate specific features in 
the baryonic matter (Malaney & Butler 1989). Yokoyama 
& Sato (1991), Dolgov & Silk (1993), Polarsky & Starobin- 
sky (1994), Novikov, Schmalzing & Mukhanov (2000) have 
shown a few different ways for the generation of the isocur¬ 
vature perturbation in the framework of the inflation theory. 

In connection with the above-mentioned problem of 
the baryonic fraction determination from the cosmological 
nucleosynthesis and the CMB anisotropy data we dedicate 
our paper to the re-examination of the models with non¬ 
linear sub-horizon-scale (at the epoch of the hydrogen re¬ 
combination) baryonic clouds with 5 b 3> 1. Such inhomo¬ 
geneities do not manifest in the CMB anisotropy( because 
of the extremely small scales) but manifest themselves by 
the transformation of the ionization history of the primeval 
hydrogen-helium plasma at redshift z ~ 10 3 . All these fac¬ 
tors should be taken into account in the reconstruction of 
the ionization history of the Universe especially at the pe¬ 
riod of the cosmological hydrogen recombination . The rea¬ 
son for importance of the possible very small scale entropy 
perturbations at the epoch of recombination is connected 
with the very small mean free path of the Ly-a photons at 
2 ~ 10 3 : 1 Lq ~3x 10 10 ((l + 2)/1000)“ 5/2 (D(,h 2 /0-02)- 1 cm 
which corresponds to the baryon mass Mlu = 4-7r/3p6Z 3 a ~ 
2 x 10 _23 ((l + z)/1000) _9/2 (n6/i 2 /0.02) _2 M Q where M 0 is 
the Solar mass. High amplitude baryonic clouds with masses 
M Mlc could transform the process of recombination 
at the beginning and dissipate during recombination up to 
the crucial masses Mdiss ~ 10 s Mq (Liu et al. 2001). We 
will show that if the typical masses of the clouds M are 
M > Mdiss the hydrogen and helium recombination inside 
and outside clouds goes independently. Due to non-linear 
dependency of the electronic ionization fraction x e on the 
baryon density the hydrogen and helium inside the clouds 
recombine faster than outside them. Thus the dynamics of 
the mean value of the electronic ionization fraction x e which 
plays a crucial role in the CMB anisotropy and polarization 
formation decreases slower than, for example, in the uni¬ 
form model with the mean value of the baryonic fraction 
of the matter. We will show that in the cloudy baryonic 
plasma the kinetics of the H — He 4 recombination is closer 
to the delayed recombination model by Peebles et al. (2000) 
with concrete relation between e 0 and a parameters of their 
model and amplitudes of perturbations and the filling factor 
of the baryonic clouds. We will show how sensitive the Ci 
power spectrum of the CMB anisotropy is to the mentioned 
above parameters of the baryonic clouds. We will discuss 
possible manifestation of the small scale perturbations in 
the MAP and upcoming PLANCK observational data. 
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2 BASIC DEFINITIONS AND 
MODIFICATIONS OF THE 
HYDROGEN-HELIUM IONIZATION 
HISTORY. 

We consider a model with non-linear baryonic perturbations 
at the small scales ( M <C 10 1o Mq). In the analysis of the ki¬ 
netics of recombination we will take into account electrons, 
protons, ionized and neutral hydrogen and helium. For sim¬ 
plicity we suppose that all baryonic clouds have the same 
characteristic sizes R c i, which are much smaller than the 
size of the horizon R rec close to the period of recombination 
(z ~ 10 3 ), Rd <C Rre c- We denote pb,in, pb,out and pb.mean 
the baryon density inside the clouds, outside of them and 
the mean density at the scales much greater than R c i and 
distances between them, respectively. We have the following 
relations 

Pb,mean — Pb,inf Pb,out (1 -/), (2) 

where / is the volume fraction of the clouds. We denote 

?= PM IL _ ( 3 ) 

Pb,out 

We can write down the following relations between 
mean value of the baryon density and inner and outer values 


ipb,mean 

Pb ’ m 1 + /(£ - 1) ’ 

(4) 

and 


Pb,mean 

Pb,out - J + /( £ _ 1} • 

(5) 


As mentioned in Introduction the presence of the bary¬ 
onic clouds in the primordial hydrogen-helium plasma at 
redshift z ~ 10 3 changes the dynamics of the recombination 
due to non-linear dependence x e on the baryon density. Be¬ 
low we will describe the kinetics of the recombination in a 
cloudy baryonic fraction of the Universe taking into account 
that the diffusion damping is not important for the clouds 
with M > Mj. As it was shown by Liu et al. (2001), during 
the period of recombination diffusion of baryons from inner 
to outer regions of the clouds can suppress any small scale 
irregularities inside the clumps. The natural length of this 
process is close to the Jeans length Rj ~ c s ri rec where c s 
is the baryonic speed of sound and r/ rec is the correspond¬ 
ing time when the plasma became transparent for the CMB 
radiation. Note that our aim is to predict some possible ob¬ 
servational features in the CMB anisotropy and polarization 
power spectrum connected with possible non-uniform distri¬ 
bution of baryons on very small scales: much smaller than 
Rrec, but greater than the Jeans scale Rj. That means that 
for all adiabatic perturbations at the scales M 10 13 Mq 
( which are the source of the Doppler peaks in the CMB 
anisotropy and polarization power spectra ) the evolution 
during the period of recombination depends not on the ion¬ 
ization fraction inside or outside the clouds but rather on the 
mean value of ionization fraction over the scales of adiabatic 
perturbations. As we mentioned above this mean ionization 
fraction does not correspond to the ionization fraction for 
the mean value of the baryonic density due to non-linear 
effects. For our analysis we use two basic software pack¬ 
ages, RECFAST (Seager et al. 2000) and CMBFAST (Sel- 
jak & Zaldarriaga 1996) with modification for the cloudy 
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baryonic model. Let us start from CMBFAST modification 
because the calculation of the CMB anisotropy and polar¬ 
ization power spectra depends on the number density of free 
electrons. 

For the cloudy baryonic model we introduce the mean 
value of the electron density 


(n e ) = n e ,inf + n B , ou t(l - /), 


( 6 ) 


where / is the fraction of volume with clouds from Eq.(2). 
For the baryonic clouds with scales R > Rj we can neglect 
diffusion of baryons and Ly — a photons and describe the 
recombination process inside and outside the clouds sepa¬ 
rately. In such a case we can write down 


( 1 1 He,in \ 

Tie,in — %e,in I 1 “ ) Tib,in’) 

f 1 YHe,out \ ^ 

Tie,out — %e,out l J- ~ J Tib,out) 


(7) 


where x e ,i n , x e ,out, Y H e,in and Y He ,out are the ionization 
fraction and helium mass fractions for inner and outer re¬ 
gions. Note that by definition x e = n e /riH, where nrr is 
the number density of neutral and ionized hydrogen. Let us 
introduce the mean value of the ionization fraction (x e ), 


{x e ) 


(ne) ( (Y He )\ 

(n b ) V 2 J 


then 


(x e ) = x e,in G in T Xe,outGouti 


( 8 ) 

(9) 


where 



Gin — 

a 

( 1 — lire, in/2 \ 

i + m - 1) 

^ 1 - (Yrre)/ 2 ) 

Gout — 

!-/ / 

f 1 — lire, oiit/2 \ 

l + /(£-D \ 

v 1 - (Y H e)/2 J 


and (Y/fe) denotes the mean mass fraction of helium. 

The second remark is related with the modification of 
the CMBFAST code, particularly with the characteristic 
time of friction td between electron and radiation fluids be¬ 
fore and at the period of the recombination. Let us consider 
the hydrodynamic equations for baryon-electron fluid and 
radiation using the Newtonian approximation . Following 
Liu et al. (2001), we have 


Ob + 3 H pb + S7i{pbVi,b) — 0, 

V iib + HVi,b + Vj,b(V jVi,b) + V i Pi, j pb + Vf'F 
4 p 7 on e <TT 


3 Oh 


-(Vi,-y-Vi, b ), 


(11) 


where H = o' /a, a is the scale factor of the Universe, 'F is 
the gravitational potential, dot denotes the derivative with 
respect to time r/ (dr/ = dt/a , the speed of light c = 1), 
p 7 is the density of the CMB, n e is the local concentra¬ 
tion of electrons, Vi , 7 is hydrodynamic velocity of radiation 
( dipole moment) and or is the Thompson cross-section. In 
our cloudy baryonic model we can define the large scale adi¬ 
abatic tail of the perturbations of the matter and the small 
scale CIM tail as follows, 


pb = pb,s(r, t) + pb,i(r, t), 

Vi, b = + 

® + ( 12 ) 
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where 


{pi) 

= (pb, s {r,t)) + p b ,i(r,t), 

(p b ,s(r,t)) 

— pb,mean(t ), 

Pb,i(r,t) 

= pb,mean(t) (1 + 5(r, t)) , 

(Vi,b) 

= (^b ) (u*)) + y i!b(u t )> 


= Hn, 


= (^ s (t)) + ^(r,i), 


and vY(r,t) and S^i{r,t) are the linear adiabatic 

perturbation of the hydrodynamic quantities. By the def¬ 
inition (...) means average over the scales |r| 2> R c i but 
R cl \k\ <C 1 for all k Fourier harmonics of the adiabatic per¬ 
turbations. Using Eq. © and © we can describe the evo¬ 
lution of the velocity perturbations for the adiabatic modes 

dV l 

+ HVl h + S(ViP to t/ptot) + 

or] 


= ( Apia 3 n p e b aT )(vL-vl b ), 


(14) 


where Ptot and ptot are the pressure and the density of the 
baryon-photon fluid and <5(V iPtot/ptot) means a linear part 
of the perturbations. As one can see from Eq. (jlj) in the 
cloudy baryonic matter the characteristic time of baryon- 
radiation friction is 


/4 p-yaneCr \ 4p 7 a(Xe)<7T 


3 pb 


'-) = 


2>m„ 


(15) 


where m p is the proton mass. This time td depends on the 
local ratio (x%) = (n e (f)/riH(r))■ The ionization fraction 
(xt) could be written down in the following form, 


(xl) 


Xe,inQin ~\~ Xe,outQout 5 


(16) 


where 


9 in 


f 1 — YHe,in/2\ 

\ 1 - {Y He )/2 ) : 


Qout 


(1 -/) 


/ 1 — Fife,out/2\ 

^ 1 - (Tife)/2 ) ■ 


(17) 


Thus Eq. ( 0 ) and Eq. ( 0 ) generalize the standard ho¬ 
mogeneous model of the CMB anisotropy and polarization 
formation on the cloudy baryonic model of the Universe. 
Note the difference between Eq. (|l4j) and the standard ki¬ 
netic equation for the CMB anisotropy formation. Whereas 
Eq. (jl4f) depends on (xi)), the standard equation depends on 
(x e ),(see Eq.(8) and (9)). The corresponding modification 
of the CMBFAST code leads to modification of RECFAST 
programme which calculates the kinetics of recombination 
inside and outside the small-scale baryonic clouds as a func¬ 
tion of cosmological parameter (fit,) = pb,mean I per-, density 
contrast £ and volume fraction /. In Fig.l and Fig.2 we plot 
the functions x e ,i„, x e , 0 ut, (x e ), and (x*) for different mean 
values of the baryonic density (fib) and h = 0.65, £ = 11, 
and / = 0.1. As one can see from these figures the frac¬ 
tions of ionization (x e ) and (x%) have different shapes and 
different asymptotic at low redshifts. 

In the next section we describe the results of the cor¬ 
responding numerical computations. However, some prelim¬ 
inary discussions could be very useful for understanding of 
the most important peculiarities connected with two char¬ 
acteristic functions of ionization (x e ) and (xl). Firstly we 



z+1 


Figure 1. The fraction of ionization ( x e ) as the function of z for 
different mean values of the baryonic density (fib) and £ = 11, 
/ = 0.1. From bottom to top three solid lines correspond to the 
fractions of ionization for inner regions of the clouds at ( flb)h 2 = 
0.03, ( flb)h 2 = 0.02 and ( fl b )h 2 = 0.01 (h = 0.65). Dash lines 
correspond to the same but for the outer regions and marked 
lines correspond to ( x e ) for the same models. 



z+1 


Figure 2. The fractions of ionization ( x e ) and (x)) as the func¬ 
tion of z for J = 11, and / = 0.1. From bottom to top three solid 
lines correspond to the fractions of ionization (x e ) from Fig.l for 
(1 fl b )h 2 = 0.03, (fl b )h 2 = 0.02 and ( fl b )h 2 = 0.01 (h = 0.65). 
Dash lines from the bottom to top correspond to ionization frac¬ 
tions (a:*) at the same numeration. 


would like to point out that in the case when the volume 
fraction / is small (/ <IC 1, but £/ ~ 1!) the function (xl) 
should be very close to Xe.outPout while (x e ) practically does 
not change. From Eq. (|16j) and Eq. ([17]) we can immedi¬ 
ately find an asymptotic of the function (x e ) at low redshifts: 
(x e ) ~ (xt)G 0 ut/gout = {xt)/(l + /(£ - 1)) < (x*). That 
means that the Silk damping scale which is very sensitive 
on the parameter r/J 1 , which is proportional to (x)), should 
increase due to the increasing of the effective ionization ratio 
(x|). Thus we can conclude that in cloudy baryonic model 
the position and amplitudes of the Doppler peaks in the Ci 
power spectrum differ from the same values in the standard 
non- cloudy model with the same mean value of the baryon 
density (fib). 

Secondly, it is well known that the so-called “visibil¬ 
ity function” g(r) = r'exp(—r), where r is the Thomp- 
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son optical depth, depends on the function (x e ). Because 
(x e ) < (**), and {xl) ~ x 0 ut we can conclude that in our 
model the kinetics of recombination is similar (in general) 
to the models with additional sources of the ionization. 

One additional comment is related with the delayed re¬ 
combination model by Peebles et al. (2000) mentioned in 
Introduction. We would like to note, that the physical basis 
of their model and our cloudy baryonic model are completely 
different, but the numerical data for the corresponding func¬ 
tion Xe in the model by Peebles et al. (2000) and ( x e ) in 
our model are close to each other. If we compare the solid 
lines with the marked lines in Fig.l, we can conclude that 
change of the ionization fraction ( x e ) plays a role in the cor¬ 
responding delay of the recombination with respect to the 
recombination inside the clouds. To understand the change 
of the CMB anisotropy and polarization power spectra we 
need to compare the ionization fraction for outer regions and 
(x e ). From Fig.l it is clearly seen that in such a case we can 
introduce the term “accelerated recombination” because the 
mean ionization ( x e ) —> 0 faster then x e —> 0 for the outer 
regions at the same values of the Llbh 2 parameter. In both 
cases, following Peebles et al. (2000), we can describe the 
number of additional resonance Ly — a quanta produced at 
the epoch of the hydrogen recombination by the sources of 
ionization 


dvij-es 

dt 


eH{t)rib , 


(18) 


where e is the efficiency of the Ly — a quanta produc¬ 
tion, H{t) is the Hubble parameter and rib is the num¬ 
ber density of baryons. According to Peebles et al. (2000), 
e > 0 describes the productivity of the sources of additional 
resonance quanta, which leads to the delay of the hydro¬ 
gen recombination. Formally, if e < 0 then recombination 
goes faster and we have “accelerated recombination” regime. 
Thus, the difference between these two situations is related 
to the definition of the background state (inner or outer 
parts of the clouds). Note that for the CMB power spec¬ 
trum calculations the term “accelerated recombination” is 
preferable because it depends mainly on the characteristics 
of the outer zones. 


3 ANISOTROPY AND POLARIZATION 
POWER SPECTRA IN A CLOUDY 
BARYONIC UNIVERSE. 

For numerical calculations of the CMB anisotropy and po¬ 
larization power spectra we will use the modified CMBFAST 
code taking into account the above-mentioned peculiarities 
of the ionization history of the plasma. We take into account 
the difference of the He 4 mass fractions Yh b for inner and 
outer zones into account, using well estimated dependence of 
Yh e on (f ph 2 ) parameter (Olive, Steigman & Walker 2000). 
For illustration of the importance of the cloudy structure at 
small scales on the CMB anisotropy formation in Fig.3 we 
plot the function A 2 T(l) = l(l + l)Ci/2n (/. iK 2 ), where Ci is 
the anisotropy power spectrum, for the cosmological model 
with (Ll c dm) = 0.3, Ha — 0.65, h = 0.65, (Hb)ft 2 = 0.02 and 
A 2 T(l) for the corresponding model with the uniform bary¬ 
onic matter distribution and scale invariant power spectrum 
of the initial adiabatic perturbations. We choose / = 0.1, 
and density contrast between inner and outer zones £ = 11. 



Figure 3. The A T 2 (l) function as a function of multipole num¬ 
ber l for £ = 11 , / = 0.1 model. The dash line corresponds to 
the cosmological model for f2j, A 2 = 0.02 without clouds. Solid line 
corresponds to the model with clouds and {flb)h 2 = 0.02. 



Figure 4. The polarization A 2 T v (l) = 1(1 + l)C p (l)/2ir as a 
function of the multipole number l for £ = 11, / = 0.1 model. 
The numeration of the lines is the same as in Fig.3. 


As one can see from Fig.3, the presence of the clouds 
before and during the period of the hydrogen and he¬ 
lium recombination significantly perturbed the correspond¬ 
ing A T 2 (l) function. The same conclusion follows from Fig.4 
for the polarization of the CMB in the cloudy baryonic 
model. 

As it is seen from Fig.4 the more complicated ioniza¬ 
tion history of the plasma in the cloudy model leads to the 
decreasing of the corresponding power spectrum of the po¬ 
larization at l > 900 due to the increasing of the effective 
dissipation scale. It is worth noting that in our cloudy model 
we do not take into account the possible reionization of the 
hydrogen at low redshift (z < 20) due to the influence of 
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Figure 5. The functions D(l) for the anisotropy (solid line) and 
polarization ( dashed line) as functions of the multipole number 
l for £ = 11 , / = 0.1 model. 

the additional sources on the ionization balance. This mod¬ 
ification is the standard part of the CMB anisotropy and 
polarization spectrum calculations using CMBFAST codes. 
But it is clear that in cloudy baryonic Universe the mean 
value of the optical depth must be related to the mean ion¬ 
ization fraction (x e ) for corresponding values of the (fij,) 
parameters. 


4 CONCLUSION. 


can see from Fig.5, practically for all ranges of multipoles 
the differences between cloudy and non-cloudy models are 
observable for the PLANCK mission. Moreover, if the pa¬ 
rameter £ > 1, then we need to include the possible cloudy 
baryonic model in the schemes of the cosmological parameter 
extractions from the current and future CMB anisotropy and 
polarization data. As one can see from Eq.(4) and Eq.(5), 
if /(£ — 1) <C 1, then the difference between pb,out and 
mean baryon density is ~ /(£ — 1). According to the pre¬ 
dicted accuracy of the cosmological parameter extraction 
from the PLANCK mission, the corresponding uncertain¬ 
ties for the baryonic fraction of the matter Aft = 8Qb/^b 
must be less then a few percents. Taking conservative limit 
At ~ /(£ — 1) ~ 0.1 and £ — 1 ~ 1 we can obtain that 
the corresponding fraction / should be detectable by the 
PLANCK satellite, if / > 0.1. 

We would like to point out that our simple model of 
the baryonic clouds is based on the one possible modes of 
the isocurvature perturbations at the small scales. It would 
be interesting to investigate the more complicated models of 
the initial perturbations. This program is in progress. 
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As it was mentioned in Introduction, the baryonic fraction 
of the matter is one of the most important cosmological pa¬ 
rameters, which determines the most preferable model of the 
Universe. Two crucial parameters of the theory are now un¬ 
der discussion, i.e., the light chemical elements abundance, 
which is related directly to the SBBN predictions, and the 


Iff, parameter which rail be determined from the current and 
future CMB observations. In both cases the determination 
of the parameters depends on the hypothesis about the spa¬ 
tial distribution of the baryonic fraction of the matter and 
can be tested by the CMB experiments such as the launched 
MAP and the future PLANCK missions. 

In our paper we have shown how important it is the pos¬ 
sible baryonic inhomogeneity at small scales (M <C 10 13 Mq) 
in the CMB anisotropy and polarization formation through 
distortions of the cosmological recombination. The extremal 
case, when the density contrast £ ~ 10, shows that the 
clumps in the baryon fraction of the matter can significantly 
change the amplitudes and positions of the Doppler peaks 
in the CMB anisotropy and polarization power spectrum. 
For example in the above mentioned £ ~ 10 model the dif¬ 
ferences between Ci for the anisotropy and for the polar¬ 
ization in the models with flbh 2 = 0.02 ( no clouds) and 
(! Qb)h 2 = 0.02 (with clouds) are presented in Fig.5 in a form 
of the following the functions. For the anisotropy 


D a (l) = 2{Ca,nc{l) - Ca,oil))/{C a ,noil) + C7„, c (J)), 


(19) 


where C a ,nc{l) and C a ,c(l) denote the non-cloudy and cloudy 
model, respectively, and the index a corresponds to the 
anisotropy spectrum. The analogous definition of the D p (l) 
function is used in Fig.5 for the CMB polarization. As one 
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